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摘要 :生态 系统 碳 循 环 过 程 对 水 分 响应 的 研究 已 成 为 全 球 变化 关注 的 焦点 问题 之 一 。 基 于 长 白山 温带 针 闪 混交 林 与 千 烟 洲 亚 
热带 人 工 针 叶 林 观测 站 2003 一 2009 年 生长 季 的 碳 通 量 (NEE ) 和 气象 观测 数据 ,综合 考虑 水 分 对 光合 ,呼吸 作用 的 影响 ,构建 不 
同 的 NEE 模型 ,并 应 用 模型 数据 融合 方法 优化 模型 参数 . 送 选 最 适 模型 ,系统 分 析 了 水 分 因子 对 不 同 森林 生态 系统 碳 循环 的 影 
响 。 结 果 表 明 :(1) 优 化 后 的 模型 参数 均 能 被 NEE 实测 数据 较 好 约束 。 长 白山 生长 季 的 光合 .呼吸 参数 值 均 高 于 千 烟 洲 , 未 考 
虑 空气 饱和 水 汽 压 差 (VPD) 的 模型 高 估 了 千 烟 洲 温度 敏感 性 参数 (Qio ) 值 .低估 了 千 烟 洲 基础 呼吸 速率 参数 (8) 值 ;(2) 仅 考 
虑 VPD 对 光合 作用 影响 的 模型 是 长 白山 生长 季 碳 通 量 模拟 的 最 优 模 型 ,但 模拟 精度 提高 不 显著 。 不 同 模型 间 碳 通 量 组 分 模拟 
结果 差异 较 小 ;(3) 考虑 VPD 和 土壤 含水 量 对 光合 .呼吸 作用 共同 影响 的 模型 是 千 烟 洲 生长 季 碳 通 量 模拟 的 最 优 模型 ,并且 显 
著 提 高 了 模拟 精度 。 未 考虑 水 分 的 模型 在 生长 季 高 估 了 总 生态 系统 生产 力 (GEP) 总 量 2.0%(21.85 g C/m ) ,同时 更 大 幅度 地 
高 估 了 生态 系统 呼吸 ( RE) 总 量 4.4%(38.02 g C]m  ) ,从 而 导致 NEE 总 量 低估 于 实测 值 7.8%(18.55 g CA/m?)。 
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Abstract: Moisture effect on the carbon balance of terrestrial ecosystems is a key issue in global change research. It is 
crucial to accurately analyze the response of terrestrial ecosystem carbon cycle to moisture. However, the carbon flux models 
responding to environmental factors rarely consider the moisture effects on photosynthesis and respiration simultaneously; 


meanwhile there are still large uncertainties in model structures and parameters. Thus, this study was designed to (1) 
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choose the optimal carbon flux model with accurate parameters for different ecosystems through model-data fusion approach, 
reducing the uncertainties of modeled results; (2 ) systematically analyze the influence of water factors on carbon flux 
simulation, including gross ecosystem productivity ( CEP ) ，ecosystem respiration ( RE) and net ecosystem exchange 
(NEE ). To consider the effects of moisture on both GEP and RE, we developed four different NEE models. Then, based on 
carbon flux and meteorological data during growing season from 2003 to 2009 in Changbaishan temperate mixed forest 
(CBS) and Qianyanzhou subtropical coniferous plantation (QYZ ) Markov Chain Monte Carlo was employed to estimate 
model parameters, and Bayesian Information Criterion was applied to choose the optimal model for two forest ecosystems. 
The results showed that (1) the posterior values of model parameters were normally distributed, indicating that the 
parameters were well constrained by NEE. Photosynthetic and respiratory parameter values of CBS were higher than those of 
QYZ during the growing season. The model without vapor pressure deficit ( VPD) overestimated the value of temperature 
sensitivity ( O,) and underestimated the value of basal respiration rate ( B,) in QYZ; (2) the model considering VPD only 
was the optimal model for CBS, but its performance was not improved much. The modeled flux components were similar 
among the four models; (3) the model considering both VPD and soil water content (S$,) was the optimal model for QYZ ， 
and its performance was improved significantly. The model ignored water factors overestimated 2% (21.85 g C/m ) of the 
total GEP, and 4.4% (38.02 g C/m ) of the total RE, and therefore, underestimated 7.8% (18.55 g C/m’) of the total 


measured NEE during the growing season. 


Key Words: moisture effect; carbon cycle; model-data fusion; parameter optimization; model selection 


模型 是 理解 碳 循 环 过 程 及 其 控制 机 制 的 重要 手段 ' 中 ,如何 准 确 地 模拟 、 预 测 陆 地 生态 系统 碳 循环 变化 过 
程 成 为 目前 全 球 变化 最 为 重要 的 研究 问题 之 一 ' 站 。 通 量 观测 技术 及 全 球 观测 网 络 的 发 展 ,为 陆地 生态 系统 
碳 循环 模拟 研究 提供 了 大 量 连续 观测 的 数据 支持 '* | ,光合 .呼吸 作用 对 光照 温度、 水 分 等 环境 要 素 变化 的 
响应 模型 ( 碳 循 环 过 程 环 境 响 应 模型 ) 也 随 之 得 到 了 长 足 的 发 展 ,成 为 跨 尺 度 生 物 地 球 化 学 循环 模型 的 重要 
组 成 部 分 ,在 陆地 生态 系统 碳 循环 过 程 模 拟 中 起 着 重要 的 作用 “7 。 水 分 是 影响 生态 系统 碳 循环 过 程 的 重要 
因素 ,尤其 是 那些 水 分 受 限 的 生态 系统 !5 。 它 通过 影响 微生物 活性 、 植 物 根系 生长 .营养 物质 迁移 及 呼吸 的 
温度 敏感 性 而 限制 生态 系统 呼吸 作用 !”) ,也 通过 调节 植物 气孔 导 度 和 叶片 水 势 . 影 响 光 合 系统 的 酶 活性 和 
C0, 固 定 效率 而 限制 生态 系统 光合 作用 ' 咏 ] 。 目 前 在 全 球 变化 背景 下 ,降水 变 率 的 增 大 导致 干旱 影响 区 域 和 强 
度 的 增 大 "由 ,生态 系统 碳 循 环 过 程 对 水 分 响应 的 研究 越 来 越 得 到 关注 和 重视 "1 。 其 中 Richarson 等 [| 、Liu 
等 中 构建 呼吸 作用 对 温度 .水 分 变化 的 响应 模型 ,系统 分 析 了 水 分 对 生态 系统 呼吸 的 影响 ;Lassop 等 1、 
Wagle 等 09 基于 多 个 通 量 站 数据 构建 光合 作用 对 光照 .空气 湿度 变化 的 响应 模型 ,探讨 了 水 分 在 不 同 植被 光 
合作 用 中 的 重要 性 。 但 水 分 变化 同时 影响 生态 系统 的 光合 .呼吸 作用 ,而 现 有 碳 循环 过 程 环境 响应 模型 较 少 
考虑 水 分 因子 对 光合 .呼吸 的 共同 影响 ,无 法 准确 分 析 水 分 对 生态 系统 碳 循环 过 程 的 影响 。 

碳 循环 过 程 对 水 分 响应 形式 的 不 同 使 模型 结构 存在 差异 性 ,同时 模型 参数 也 难以 准确 估算 ,导致 模拟 结 
果 存 在 较 大 的 不 确定 性 。 而 近年 来 发 展 的 基于 贝 叶 斯 理论 的 模型 数据 融合 方法 (Model-Data Fusion ， 
MDF) 3 能够 充分 利用 观测 数据 和 参数 先 验 信息 ,优化 模型 参数 .评估 模型 结构 ,为 准确 分 析 水 分 对 生态 系统 
碳 循环 影响 提供 了 有 效 的 技术 手段 !'”) 。 目 前 MDF 在 碳 通 量 模拟 中 已 得 到 广泛 应 用 ,但 主要 集中 在 通过 马尔 
科 夫 链 -蒙特 卡 罗 (MCMC ,Markov Chain Monte Carlo) 等 优化 方法 降低 模型 参数 不 确定 性 方面 。 综 合 考虑 
模型 参数 及 模型 结构 不 确定 性 的 研究 较 少 "1 ,其 中 Sacks 等 2 、Zobitz 等 52 采用 贝 叶 斯 信息 准则 法 (BIC ) 定 
量 比较 了 不 同 呼吸 .光合 模型 的 优势 和 不 足 , 揭 示 了 MDF 应 用 在 模型 结构 评价 和 选择 中 的 优越 性 。 

于 烟 洲 亚热带 人 工 针 叶 林 (QYZ) 和 长 白山 温带 针 阔 混交 林 (CBS ) 分 别 代表 了 我 国 东 部 南北 森林 样 带 上 
的 两 个 典型 森林 生态 系统 类 型 "3 ,并 在 生长 期 水 分 条 件 差异 显著 ,这 为 研究 水 分 对 于 不 同 生态 系统 碳 循环 过 
程 的 影响 提供 了 天 然 的 实验 条 件 。 其 中 QYZ 雨 热 不 同 季 造成 了 季节 性 干旱 ,对 生态 系统 碳 循 环 过 程 存在 胁 
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迫 现象 ,而 CBS 在 生长 季节 水 分 相对 充沛 , 较 少 受到 干旱 的 影响 。 因 此 本 文 以 QYZ 和 CBS 为 例 ,综合 
考虑 水 分 对 光合 .呼吸 作用 的 共同 影响 ,构建 不 同 的 碳 循环 过 程 环境 响应 模型 ,集成 MCMC 与 BIC 方法 优化 
模型 参数 .评价 模型 结构 , 遂 选 两 个 生态 系统 的 最 适 模 型 ,降低 碳 通 量 模拟 结果 的 不 确定 性 ;并 在 此 基础 上 探 
讨 水 分 因子 对 QYZ .CBS 光合 .呼吸 作用 的 影响 。 


1 数据 与 方法 


1.1 站 点 概况 和 研究 数据 

千 烟 洲 (QYZ) 观测 站 位 于 江西 省 红壤 丘陵 农业 综合 开发 实验 站 内 ,多 为 树龄 约 30a 的 人 工 针 叶 林 , 属于 
典型 亚热带 季风 气候 ;长 白山 (CBS) 观 测 站 位 于 长 白山 自然 保护 区 红 松 针 阔 混交 林内 ,多 为 树龄 约 200a 的 成 
熟 原始 林 , 属于 典型 的 温带 大 陆 性 季风 和 气候。 南北 两 个 生态 系统 水 热 状 况 差 异 明显 。 中 国生 态 系统 研究 网 络 
(CERN) 和 中 国 陆 地 生态 系统 通 量 研究 网 络 ( ChinaFLUX) 对 两 个 站 点 都 进行 了 长 期 的 生态 监测 和 碳 通 量 连 
续 观 测 ,站 点 的 基本 特征 见 表 1 ,其它 详细 信息 可 以 参考 张 雷 明 等 .3 刘 允 芬 等 .| 。 


表 1 研究 站 点 的 基本 特征 


Table 1 Basic characteristics of research sites 


站 点 Site 千 烟 洲 亚热带 人 工 针 叶 林 ( QYZ) 长 白山 温带 针 阔 混交 林 ( CBS) 
Qianyanzhou subtropical coniferous plantation Changbaishan temperate mixed forest 

地 理 位 置 Geographical location 26°44'N,115°04'E 41°29'N,128°05'E 

海拔 Elevation /m 100 736 

年 平均 温度 Mean annual temperature /% 17.9 4 

年 降水 Mean annual precipitation /mm 1485 695 

土壤 类 型 Soil type 红壤 山地 上 暗 柠 壤 

冠 层 高 度 Canopy height /m 11 26 


研究 数据 包括 模型 驱动 . 反 演 及 验证 数据 。 模 型 驱动 数据 为 两 个 站 点 2003 一 2009 年 的 30 min 气象 观测 
数据 ,主要 包括 气温 .土壤 温度 土壤 含水 量 、 光 合 有 效 辐 射 .空气 饱和 水 汽 压 差 , 均 由 常规 自动 气象 观测 系统 
测定 。 反 演 和 验证 数据 均 为 两 个 站 点 经 过 质量 控制 后 的 30 min 涡 度 相关 碳 通 量 观测 数据 ,具体 处 理 流程 参 
见 李 春 等 :51 。 其 中 ,2003 一 2006 年 生长 季 (5 一 9 月 ) 的 通 量 数据 用 于 模型 参数 优化 ,2007 一 2009 年 生长 季 
(5 一 9 月 ) 的 通 量 数据 用 于 模型 验证 。 
1.2 NEF 模型 组 合 
净 生 态 系统 碳 交 换 量 ( NEE) 是 总 生态 系统 生产 力 ( GEP ) 与 生态 系统 呼吸 ( RE) 两 大 通 量 平衡 的 结果 。 
本 研究 在 温度 响应 的 呼吸 模型 与 光 响 应 的 光合 模型 的 基础 上 ,综合 考虑 土壤 水 分 和 空气 水 分 因子 对 GEP .RE 
的 影响 ,构建 了 4 组 NEE 模型 表述 碳 循环 过 程 对 水 分 响应 的 不 同形 式 。 表 2 描述 了 碳 通 量 模 型 的 组 合 方 式 。 
Lloyd&Tayor 模型 ”1( 1) 作 为 生态 系统 呼吸 的 温度 响应 模型 得 到 广泛 认可 和 应 用 , 相 比 Van' Hoff 方程 和 
Arrhenius 方程 , 它 更 好 地 描述 了 生态 系统 呼吸 的 变异 性 ,在 低温 和 高 温 条 件 下 对 温度 变化 都 有 较 强 的 响应 能 
力 ' 引 。0,, 模 型 (2) 采 用 温度 和 土壤 含水 量 作 为 模型 驱动 变量 来 描述 温度 和 水 分 对 生态 系统 呼吸 的 协同 作 
用 ,并 假设 生态 系统 呼吸 的 温度 敏感 性 因子 016 与 土壤 水 分 状况 密切 相关 。 相 比 于 其 他 考虑 温度 、 水 分 响 
应 的 呼吸 模型 (如 连 乘 模型 ) ,01, 模 型 在 干旱 土壤 条 件 下 对 温度 和 水 分 有 着 更 强 的 敏感 性 '*]。 
RE =Bi x QW Or (1) 
RE =B, x Ol I O00 =atb xSy, (2) 
式 中 ,RE 是 生态 系统 呼吸 (umol ms  ) ;Bi 是 参考 温度 下 的 生态 系统 基础 呼吸 (jymol ms ) ;0 是 生态 系 
统 呼吸 的 温度 敏感 性 因子 ,代表 温度 升 高 10% 生态 系统 呼吸 的 相对 增长 量 ;7 是 5cm 土壤 温度 (% ) ;7 是 参 
考 温度 (% ) ,在 此 取 10Y ;5% 是 5cm 土壤 含水 量 (m m”) ;a,b 是 试验 参数 , 当 6 为 正 时 ,表示 生态 系统 呼吸 
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对 温度 的 敏感 性 随 水 分 的 增加 而 增加 。 

生态 系统 光合 作用 强度 与 光合 有 效 辐射 旦 显著 相关 关系 ,通常 符合 直角 双 曲 线 方程 ,可 以 用 Michaelis- 
Menten 模型 (3) 来 表达 !594 。 当 考虑 GEP 对 光合 有 效 辐射 和 饱和 水 汽 压 差 的 协同 响应 时 ,可 在 Michaelis- 
Menten 模型 的 基础 上 将 4,., 蔡 换 成 一 个 考虑 VPD 的 指数 递减 函数 (4)' ,使 4 在 VPD 值 较 大 时 更 符合 生 
态 学 意义 , 即 空气 越 干 燥 , 生 态 系统 对 二 氧化 碳 吸 收 能 力 越 弱 。 


7 LUE 
EP = 一 4 
G "JT. LUE + A (3) 
4 exp( ~ k(VPD - VPD,) ) ,VPD > VPD, 
4 = a 
四 A ,VPD < VPD, ( ) 


式 中 ,GEP 是 总 生态 系统 生产 力 ( pmol ms ') ;LUE 是 初始 光 能 利用 效率 ( hmol C0,/pmol light) ; 4,, 为 最 
大 光合 速率 即 光 饱 和 时 生态 系统 CO, 同化 能 力 (pmol CO, m-? s-) ; I 为 光合 有 效 辐 射 即 入 射 到 植被 上 方 的 
光量 子 通 量 密度 (jmol light ms ) ; VPD 为 冠 层 高 度 的 空气 饱和 水 汽 压 差 (hPa) ,VPD, 设 定 为 10 hPa;4,,w 
是 不 受 VPD 约束 下 的 最 大 光合 速率 ;k 是 生态 系统 碳 吸收 对 VPD 响应 的 最 大 系数 。 
NEE=CGEP+RE (5) 
表 2 不 同 光合 呼吸 模型 组 合 


Table 2 Four combinations of photosynthetic models and respiration models 


模型 组 合 ( NEE) 呼吸 模型 ( RE) 光合 模型 (GEP) 参数 

Model combination Respiration model Photosynthetic model Parameters 

模型 1 Model 1 Lloyd&Tayor Michaelis-Menten Br Qo A LUE 

模型 2 Model 2 Lloyd&Tayor Michaelis-Menten -VPD BR OO 4 LUE 大 
模型 3 Model 3 Oo-Sw Michaelis-Menten BrQio a b A LUE 
模型 4 Model 4 Qi0-Sw Michaelis-Menten -VPD BrQio a b A LUE 大 


1.3 ”参数 估计 方法 

本 文采 用 的 马尔 可 夫 链 -蒙特 卡 罗 方 法 是 以 贝 叶 斯 定理 为 基础 ,利用 计算 机 进行 模拟 的 蒙特 卡 罗 方 法 。 
所 它 通过 构造 马尔 可 夫 链 来 推断 参数 的 后 验 信息 ,不断 迭 代 搜 索 到 全 局 最 优 解 。 相 比 最 小 二 乘法 . 极 大 似 然 法 
- 亿 等 传统 模型 参数 估计 方法 ,MCMC 方法 能 解决 包含 大 量 参数 的 复杂 非 线性 问题 ,并 结合 了 生态 模型 参数 先 验 
知识 和 观测 样本 的 信息 ,可 以 更 加 准确 地 估计 参数 。 同 时 它 还 能 提供 参数 后 验 分 布 和 置信 区 间 ,便于 进行 不 
确定 性 分 析 。Metropolis-Hastings 算法 5 是 最 常用 的 MCMC 抽样 方法 ,其 主要 步骤 为 ; 

1) 在 参数 的 先 验 范围 内 取 参 数 初 值 9,( 可 随机 也 可 指定 ) ,将 迭代 次 数 1 设 为 0 

2) 从 建议 分 布 中 产生 参数 建议 值 7; 

3) 产 生 符合 (0,1) 均 匀 分 布 的 随机 数 Ui 

4) 计算 接受 概率 a( 9,,Y) 


™(Y) (001 7) | (6) 


"7(0,) g(Y| 9,) 

如 果 a( 9,,7) 三 UV, 则 接受 参数 建议 值 , 即 9,,, =Y; 否 则 拒绝 参数 建议 值 , 即 0,,, =0,。 

5) 令 1=t+1, 重 复 2 一 5 步 , 直 到 获得 足够 多 的 样本 。 

根据 贝 叶 斯 理论 ,步骤 4 中 参数 集 后 验 概率 密度 的 求解 可 以 转化 为 似 然 函 数 的 求解 ,模型 计算 中 采用 形 
式 简单 的 对 数 似 然 函 数 , 当 参数 迭代 结束 时 便 能 获取 对 数 似 然 函 数值 的 全 局 最 优 解 。 本 文 假设 观测 数据 随机 
误差 服从 正 态 分 布 , 似 然 函 数 和 其 对 数值 分 别 如 下 : 


a(0,,Y) = al 


I 1 一 《和 一 人 了 2 
了 = e 202 (7) 
i=1 V2TO 
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LL=-N:In(V/27)-N:In(o)- 2 
式 中 ,NN 为 观测 数据 个 数 ,x, 和 ,分别 为 第 i 个 观测 值 和 模拟 值 ,o 为 观测 数据 的 标准 差 。 
1.4 模型 评价 方法 

传统 的 模型 评价 方法 包括 比较 模拟 值 与 实测 值 的 似 然 函数 值 .决定 系数 R 、 均 方 根 误差 (RMSE) 以 及 平 
均 绝 对 偏差 (MAE ) ,但 其 忽略 了 不 同 结构 的 模型 引起 的 参数 组 成 和 数目 的 差异 。 本 文采 用 的 贝 叶 斯 信息 准 
则 法 不 仅 考 虑 了 模型 数据 匹配 情况 ,还 考虑 了 参数 的 个 数 以 及 参与 拟 合 的 数据 点 个 数 , 从 而 避免 了 模型 结构 
的 宛 余 以 及 数据 的 过 度 拟 合 ,能 够 从 不 同 的 模型 中 北 选 出 结构 精简 .模拟 效果 最 佳 的 模型 ”| 。 

BIC =- 2LL + Kln(n) (9) 

式 中 ,K 是 优化 参数 的 个 数 ,n 为 优化 过 程 中 使 用 的 数据 点 个 数 ,BIC 值 越 小 说 明 模 型 越 佳 。 


2 结果 与 讨论 


(8) 


2.1 参数 优化 结果 

基于 千 烟 洲 .长 白山 两 个 森林 生态 系统 2003 一 2006 年 气象 和 碳 通 量 观测 数据 ,通过 MCMC 算法 对 4 组 模 
型 进行 参数 优化 ,获取 参数 后 验 分 布 及 最 优 参 数 集 。 图 1 和 图 2 中 4 组 模型 参数 的 后 验 分 布 均 为 正 态 分 布 ， 
且 各 参数 的 标准 差 均 较 小 ( 表 3) :呼吸 参数 标准 差 均 控制 在 参数 均值 的 10% 以 内 ,光合 参数 标准 差 均 控制 在 
参数 均值 的 1% 左 右 ,表明 优化 后 的 参数 均 可 以 被 NEE 实测 数据 良好 约束 。 


表 3 参数 优化 结果 


Table 3 Parameter optimization results 


,。 。” 先 验 值 0 风光 ,。 先 验 值 
呼吸 参数 和 模型 后 验 值 后 验 值 光合 参数 Be 模型 后 验 值 后 验 值 
Parameter value Model Posterior Posterior Parameter I Model Posterior Posterior 

value value value value 
Bp 1[ 0.25,10] 1 1.67+0.04 4.92+0.20 4 20[5,50] 1 28.96+0.37 31.58+0.25 
2 3.14+0.05 4.67+0.22 2 
3 1.84+0.05 4.96+0.21 3 29.13+0.30 31.40+0.27 
4 3.36+0.10 4.74+0.14 4 
On 2[1,3.5] 1 1.83+0.04 1.95+0.13 > 30[ 20,60] 1 
2 1.34+0.01 2.10+0.03 2 35.16+0.42 37.27+0.33 
3 3 
4 4 35.28+0.50 37.09+0.35 
a 3[1,10] 1 大 0.01[0,0.5] 1 
2 2 0.02+0.00 0.05+0.00 
3 2.22+0.03 2.41+0.18 3 
4 1.44+0.03 2.31+0.13 4 0.03+0.00 0.05+0.00 
b -4[ -10,10] 1 LUE 00. 1 0.06+0.00 0.09+0.00 
[0.001,0.1] 
2 2 0.05+0.00 0.08+0.00 
3 —1.82+0.11 -2.10+0.19 3 0.06+0.00 0.09+0.00 
4 —1.08+0.09 -1.38+0.10 4 0.05+0.00 0.08+0.00 


对 比 生态 系统 光合 参数 ( 表 3) ,不 同 模 型 间 的 参数 变化 微小 ,而 不 同 生态 系统 的 参数 差异 明显 :尽管 QYZ 
作为 人 工 林 生态 系统 具有 较 强 的 光合 能 力 ,但 CBS 在 生长 季 的 光 能 利用 效率 LUE 和 最 大 光合 速率 4,,. 仍 高 
于 QYZ。 其 中 ,LUE 平均 约 高 0.03umol CO,; pmol ! light,4， 平 均 约 高 2.20kmol CO, ms!。 这 是 因为 CBS 
温带 落叶 植物 不 同 于 QYZ 的 亚热带 常 绿 针 叶 林 , 有 明显 的 生长 季节 和 非 生长 季节 之 分 ,因此 在 有 限 生 长 季节 
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模型 1 
800 9000 9000 9000 
00 6000 6000 6000 
400 
200 3000 3000 3000 
0 0 0 0 
1.62 1.68 1.74 1.6 1.8 2.0 28.0 28.8 29.6 0.060 0.063 
BR Qno Amax LUE 
模型 2 
9000 8000 8000 8000 9000 
ebbo 6000 6000 6000 G0 
4000 4000 4000 
3000 ooo 2000 2000 3000 
位 0 0 0 0 0 
S 3.0 3.1 3.2 3.3 1.32 1.36 34.4 35.2 36.0 0.052 0.053 0.054 0.026 0.028 
六 Br Qio A LUE K 
[so 
Ee 模型 3 
芭 8000 
人 10000 8000 8000 
8000 0 6000 6000 
4000 6000 4000 4000 4000 
4000 
2000 2000 2000 2000 2000 
0 0 0 0 0 
1.76 1.84 1.92 2.15 2.20 2.25 2.30 -2.0 -1.8 -1.6 28.8 29.6 0.058 0.060 0.062 
BR a b 人 LUE 
模型 4 
8000 8000 8000 sou 2000 
en 6000 6000 e000 6000 en 
4000 4000 4000 4000 4000 
2000 2000 2000 2000 2000 3000 
0 0 0 0 0 0 
3.2 3.4 3.6 1.40 1.45 1.50 -12 -1.0 34.4 35.2 36.0 0.050 0.052 0.054 0.024 0.027 
BR a b A LUE K 
图 1 四 组 模型 参数 后 验 分 布 图 ( 千 烟 洲 和 森林 生 态 系统 ) 


Fig.1 Posterior distribution of parameters in four models at QYZ 


内 具有 高 效 利 用 有 效 水 热 资源 的 植被 生长 特性 衬 ”1 。 而 QYZ 在 夏季 受 副 热带 高 气压 带 影响 , 易 形 成 高 温和 
干旱 天 气 限制 光合 能 力 '”] 。 

对 比 生态 系统 呼吸 参数 ( 表 3) ,不 同 模型 间 的 的 参数 值 变化 较 大 ,不 同 生态 系统 间 的 参数 差异 也 明显 。 
模型 间 的 参数 变化 主要 表现 在 QYZ 的 Bi 和 01 ,CBS 的 有 和 0 变化 较 小 ,反映 出 水 分 对 QYZ 呼吸 作用 的 限 
制 。QYZ 生态 系统 的 B, 和 0 参数 集散 点 (图 3) 表 明 , 四 组 模型 B, 和 0 之 间 均 存在 很 强 的 负 相关 性 (r= 
-0.89 一 -0.97) 。 当 模型 中 考虑 水 分 因子 后 , Oo 逐步 降低 。 这 是 因为 呼吸 的 温度 敏感 性 与 水 分 呈正 相 
关 !8 ?1 ,QYZ 夏季 存在 干旱 胁迫 故 0,。 值 应 处 在 较 低 的 水 平 ,表明 考虑 水 分 后 的 模型 可 修正 Oo 的 高 估 。 
同时 ,图 3 中 考虑 VPD 因子 的 模型 2 和 4 散 点 分 布 相近 ,Bi 在 3.02 一 3.63 之 间 ,Qw 在 1.20 一 1.37 之 间 ,符合 
参考 前 人 文献 “3, 中 QYZ 呼吸 参数 值 范围 ,表明 在 QYZ 考虑 VPD 对 光合 作用 影响 的 模型 能 够 获取 
合理 的 .具有 物理 特征 意义 的 参数 ;而 未 考虑 VPD 因子 的 模型 1 和 3 散 点 分 布 相 近 : 有 范围 在 1.58 一 1.95 之 
间 ,低估 了 45%;O, 范 围 在 1.77 一 2.15 之 间 ,高 估 了 37% ,该 规律 与 Lasslop 等 !51 的 研究 结果 一 致 。 因 本 研究 
采用 NEE 数据 联合 反 演 光合 .呼吸 参数 ,两 者 存在 联动 关系 , 故 光 合 模型 中 VPD 因子 的 修订 会 影响 呼吸 参数 
的 反 演 。 呼 吸 参 数 的 生态 系统 差异 表现 为 :4 个 模型 中 CBS 生长 季 的 基础 呼吸 速率 Bi ,温度 敏感 性 参数 0 
均 显 著 高 于 QYZ。 其 中 ,Bi 平均 约 高 2.32pmol ms … ,Onu 平 均 约 高 0.44, 符 合 温带 森林 呼吸 参数 高 于 亚热带 
森林 的 规律 。 这 主要 是 因为 QYZ 的 纬度 低 于 CBS ,温度 高 于 CBS ,而 Bi 和 0Q16 均 随 着 纬度 降低 、 温 度 升 高 而 降 
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模型 1 
9000 
6000 6000 6000 
3000 3000 3000 
0 0 
4.5 4.8 5.1 5.4 0 618 2.0 3 31.0 31.5 32.0 
BR Qno Amax 
模型 2 
6000 6000 6000 
3000 3000 3000 
by 0 0 
虽 44 48 5.2 0 04 2.08 2.12 36.5 37.0 37.5 38.0 
E} 
加 BR Qno Amaxo 
区 
帮 模型 3 
紧 6000 6000 6000 
4000 4000 
3000 2000 2000 
0 0 
4.5 4.8 5.1 5.4 1 24 2.7 -2.8 -2.1 -1.4 
BR b 
模型 4 
6000 8g000 6000 
3000 4000 3000 
0 0 
44 46 4.8 5.0 20 22 24 2.6 -2.0 -1.5-1.0 
BR b 


0.090 0.095 0.100 
LUE 


8000 6000 
4000 3000 
0 


0 
0.076 0.080 0.084 0.088 0.044 0.048 0.052 


LUE K 
9000 9000 
6000 6000 
3000 3000 
0 0 
30.8 31.5 32.2 0.090 0.096 
Amax LUE 
9000 9000 
6000 
6000 6000 
3000 3000 3000 
0 0 0 
36.637.2 37.8 0.076 0.080 0.084 0.040 0.044 0.048 0.052 
Amaxo LUE K 


图 2 四 组 模型 参数 后 验 分 布 (长 白山 森林 生态 系统 ) 


Fig.2 Posterior distribution of parameters in four models at CBS 


低 '3*3”| 。 同 时 ,CBS 红 松林 林 龄 大 于 QYZ 人 工 林 , 凋 
落 物 多 .土壤 有 机 质 售 量 高 .3 故 呼吸 底 物 相 对 较 
多 ,从 而 呼吸 参数 值 较 高 。 
2.2 最 优 模 型 选择 
基于 千 烟 洲 和 长 白山 2007 一 2009 年 5 一 9 月 的 气 
象 观测 数据 ,利用 参数 优化 后 的 4 组 模型 对 两 个 森林 生 
态 系统 的 碳 通 量 分 别 进行 模拟 ,通过 BIC 方法 对 考虑 不 
同 水 分 的 模型 模拟 性 能 进行 了 定量 对 比 ( 表 4)。 在 
CBS, 仅 考虑 VPD 对 光合 作用 影响 的 模型 2 BIC 值 最 
小 ,成 为 4 组 模型 中 CBS 生长 季 碳 通 量 模拟 的 最 优 模 
型 。 相 比 QYZ 而 言 ,水 分 因子 对 CBS 模型 性 能 提升 较 
小 , 尼 仅 从 0.77 提高 至 0.83,BIC 值 降 低 了 694: 其 中 ， 
考虑 Sw 因子 使 模型 BIC 平均 仅 下 降 了 12 ,考虑 VPD 
因子 使 模型 BIC 平均 下 降 了 688 ,表明 水 分 对 光合 作用 
的 影响 更 为 明显 。 这 是 因为 CBS 属于 季风 影响 的 温带 
大 陆 性 气候 ,降水 和 温度 表现 出 一 致 的 单 峰 季 节 变 化 ， 


2.4 
” 模型 1 ” 模型 2 
^ 模型 3 。 模型 4 
S 
1.6 


1 2 3 4 
Br/(nmol m? s!) 


3 BR 和 Ci 的 最 优 参数 集散 点 图 
Fig.3 “Parameter space plots of basal respiration rate (Bn) 


Versus exponential temperature sensitivity (Oo ) 


年 均 温度 较 低 , 且 土壤 的 持 水 能 力 较 强 , 故 Sw 一 直 保持 在 较 高 的 水 平 ,对 长 白山 生态 系统 的 呼吸 作用 没有 明 


显 的 限制 作用 :3 2 。 而 由 于 CBS 在 5 月 和 9 月 的 降水 较 少 ,气温 却 保持 在 10%C 以 上 ,因此 5 和 9 月 的 VPD 
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较 高 ,6 一 8 月 的 VPD 较 低 , 形 成 了 明显 的 双 峰 季节 变化 '” , 故 VPD 对 长 白山 生长 季 的 光合 作用 有 一 定 限制 
作用 。 同 时 ,Poulter 等 ”研究 也 表明 北方 森林 生态 系统 光合 作用 对 降水 变化 更 为 敏感 ,因此 考虑 水 分 对 光合 
作用 影响 的 模型 能 够 更 为 准确 地 模拟 温带 森林 生态 系统 对 全 球 变 化 响应 。 尽 管 模型 4 的 RR 高 于 模型 2, 但 模 
型 4 中 Sw 因子 对 模拟 性 能 提升 贡献 极 小 ,反而 增加 参数 个 数 使 模型 结构 复杂 , 仅 考虑 VPD 的 模型 2 BIC 值 
最 小 ,验证 了 BIC 方法 挑选 精简 模型 的 优越 性 。 


表 4 BIC 方法 模型 选择 结果 
Table 4 Model selection results based on BIC method 


模型 水 分 因子 参数 个 数 千 烟 洲 QYZ 长 白山 CBS 

Model Moisture factor Parameters R? LL BIC R? LL BIC 
1 = 4 0.70 —18918 37876 0.77 —19685 39410 
2 VPD 5 0.80 —18557 37164 0.82 -19330 38710 
3 Sw 5 0.74 —18864 37778 0.78 —19665 39380 
4 VPD Sw 6 0.84 —18496 37052 0.83 一 19328 38716 


LL( Log Likelihood ) 为 似 然 函数 对 数值 , BIC( Bayesian Information Criterion ) 为 贝 叶 斯 信息 准则 


在 QYZ ,考虑 VPD 和 Sw 对 光合 .呼吸 作用 共同 影响 的 模型 4 BIC 值 最 小 , 且 模 型 4 中 参数 相关 性 约束 得 
更 好 (图 3) ,成 为 QYZ 生长 季 碳 通 量 模拟 的 最 适 模型 。 水 分 因子 对 QYZ 模型 模拟 性 能 有 显著 提升 ,R 从 0.70 
到 0.84 提升 了 0.14(21%) ,同时 BIC 值 降 低 了 824: 其 中 ,考虑 VPD 因子 使 模型 BIC 平均 下 降 了 719 ,考虑 Sw 
因子 使 模型 BIC 平均 下 降 了 106 ,表明 水 分 是 控制 QYZ 生长 季 碳 通 量 交 换 的 重要 环境 因子 。 这 是 因为 QYZ 
属于 典型 的 亚热带 季风 气候 ,在 夏季 受 副 热带 高 气压 带 的 控制 和 影响 ,水 热 不 同步 性 明显 , 易 形 成 伏旱 “5 。 
VPD 呈 单 峰 特征 ,在 生长 季 内 可 明显 上 升 至 2kPa, 同 时 表层 Sw 可 显著 下 降 到 0.1 mm , 故 VPD 和 Sw 对 
QYZ 生长 季 的 光合 .呼吸 作用 均 有 明显 的 限制 作用 [二 “7 。 

2.3 ”水 分 因子 对 碳 通 量 模拟 的 影响 
利用 参数 优化 后 的 4 组 光合 .呼吸 模型 对 CBS、QYZ 2003 一 2009 年 生长 季 模 型 碳 通 量 模 拟 值 进行 拆 分 ， 
并 取 其 多 年 平均 值 ( 表 5) 分 析 不 同 生 态 系 统 碳 通 量 组 分 (NEE、RE .CEP ) 对 水 分 的 响应 。 


全 表 5 不 同 模型 模拟 的 碳 通 量 组 分 生长 季 总 量 
TableSs The sum of NEE GPP and Re modeled by different models during growing season 


OO 模 理 水 分 因子 千 烟 洲 QYZ 长 白山 CBS 
Model Moisture 实测 值 模拟 值 (gC/m?) 实测 值 模拟 值 。 (g C/m?) 
factor NEE NEE GEP RE NEE NEE GEP RE 
1 全 —219.64 —1118.12 898.48 —297.56 —1158.55 860.99 
2 VPD —228.79 -1102.23 873.44 —304.23 —302.71 —1143.66 840.95 
3 Sw —225.45 —1110.08 884.63 -300.37 —1150.51 850.14 
4 VPD Sw —238.19 —235.81 一 1096.27 860.46 一 306.73 -1139.70 832.97 


在 CBS , 仅 考 虑 VPD 对 光合 作用 影响 的 最 优 模型 2 各 年 NEE 模拟 值 均 与 实测 值 最 为 接近 ,多 年 平均 的 
NEE 生长 季 总 量 模拟 值 仅 低 于 实测 值 0.5%(1.52g CAm ) 。 与 之 相 比 , 未 考虑 水 分 的 模型 1 模拟 效果 最 差 ,在 
生长 季 平 均 高 佑 了 CEP 总 量 1.3%(14.89 g Cjm2 ) ,同时 高 估 了 RE 总 量 2.4%(20.04 g C/m?) ,但 NEE 总 量 仅 
低估 于 实测 值 2.2%(6.67 g CAm?) ,表明 在 CBS 不 同 的 水 分 响应 模型 之 间 模 拟 差 异 较 小 。 

在 QYZ ,考虑 水 分 对 光合 .呼吸 作用 共同 影响 的 最 优 模型 4 各 年 NEE 模拟 值 均 与 实测 值 最 为 接近 ,多 年 
平均 的 NEE 生长 季 总 量 模拟 值 仅 低 于 实测 值 1.0%(2.38g C]m2) 。 与 之 相 比 ,未 考虑 水 分 的 模型 1 模拟 效果 
最 差 , 在 生长 季 平 均 高 估 了 GEP 总 量 2.0%(21.85 g C]m2 ) ,同时 更 大 幅度 地 高 估 了 RE 总 量 4.4%(38.02 g C/ 
m”) ,从 而 导致 整体 上 NEE 总 量 低估 于 实测 值 7.8% (18.55 g CA/m?)。 模 型 2,3 的 NEE 模拟 值 均 不 同 程度 地 
低估 于 实测 值 ,其 中 碳 通 量 组 分 对 不 同 水 分 因子 的 响应 不 同 :未 考虑 Sw 因子 的 模型 2 主要 高 佑 了 RE 总 量 
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1.5%(12.98 g C/m?) ,同时 小 幅度 高 估 了 GEP 总 量 0.5%(5.96 g C/m?)。 表 明 Sw 因子 主要 影响 了 QYZ 生态 
系统 呼吸 ,同时 也 通过 调节 冠 层 导 度 进而 影响 GEP 的 大 小 ;未 考虑 VPD 因子 的 模型 3 主要 高 估 了 GEP 总 量 
1.3%(13.81 g C]m2) ,同时 高 估 了 RE 总 量 的 2.8%(24.17 g C/m?) ,表明 VPD 因子 对 千 烟 洲 碳 通 量 组 分 模拟 
影响 更 大 。 这 是 因为 干旱 胁迫 对 千 烟 洲 生态 系 光合 作用 的 影响 较 大 “9 ,而 高 水 汽 压 差 是 光合 能 力 下 降 


的 主要 原因 中 。 同 时 由 于 千 烟 洲 人 工 林 林 龄 不 长 , 尚 处 于 较 旺 成 生长 期 ,植物 自 养 呼 吸 是 千 烟 洲 生态 系统 呼 
吸 的 主要 来 源 '”*| ,并 日 与 气温 相关 性 较 强 '”。VPD 作为 空气 温度 和 湿度 的 综合 指标 直接 影响 了 生态 系统 
呼吸 ,也 通过 影响 CEP 而 间接 影响 生态 系统 呼吸 。 

QYZ 生长 季 的 水 分 条 件 年 际 差 异 较 大 ( 丰 水 年 . 干 A 
早年 ) ,对 比 2003 一 2009 各 年 QYZ 生长 季 NEE 模拟 值 本 0 


与 实测 值 的 均 方 根 误差 (图 4) ,发 现 2003 .2007 年 考虑 
VPD 的 模型 RMSE 下 降 显著 ,其 余年 份 各 模型 间 的 模 
拟 效 果 差 异 相 对 较 小 。 而 QYZ 在 这 两 年 的 7 月份 分 别 
发 生 了 极度 干旱 和 中 度 干旱 '“ ,表明 考虑 水 分 因子 ( 尤 
其 是 VPD) 可 显著 降低 模型 在 QYZ 干旱 年 份 的 模拟 误 a 一 一 
差 。 从 QYZ 干旱 年 (2003 .2007) 生长 季 模 拟 值 与 实测 Ds 

值 的 散 点 图 (图 5) 中 发 现 ,未 考虑 VPD 的 模型 1.3 均 , , ， ， ， ， , 

有 一 个 共同 点 ,在 NEE 值 较 大 和 较 小 处 散 点 的 趋势 比 2003 2004 2005 0 2008 ”2009 

较 平缓 ,模型 模拟 存在 “高 值 上 不 去 , 低 值 下 不 来 ”的 现 

象 , 即 极 值 的 模拟 效果 要 稍 差 于 中 值 的 模拟 效果 。 尤 其 | 

当 NEE 实测 值 达 到 -20pmol m”s 时 ,模型 1、3 的 模 Fig.4 Comparison of RMSE of half-hourly modeled NEE from 
拟 值 仅 达 到 -13umol m”s 左右 ,表明 未 考虑 VPD 的 。 the four models and observed NEE at QYZ 
模型 在 干旱 时 期 低估 了 NEE 的 峰值 。 这 是 因为 干旱 期 

高 水 汽 压 差 直接 驱动 气孔 关闭 限制 生态 系统 光合 作用 ,导致 午间 生态 系统 碳 吸收 能 力 急速 下 降 王 1 ,因此 
NEE 呈现 出 非 对 称 的 日 变化 趋势 :上 午时 段 值 要 高 于 下 午时 段 ,并 在 中 午 之 前 达到 峰值 。 而 未 考虑 VPD 的 模 
型 模拟 值 则 表现 出 随 光照 变化 的 对 称 趋 势 , 从 而 在 上 午时 段 低估 NEE 下午 时段 高 估 NEE' ,并 低估 日 变化 
中 NEE 的 峰值 ,出 现 了 图 5 中 极 值 模拟 较 差 的 现象 。 考 虑 VPD 的 模型 则 可 以 准确 模拟 出 QYZ 干旱 期 NEE 
日 变化 的 非 对 称 趋势 ,因此 在 散 点 图 极 值 处 的 平缓 现象 均 有 明显 改善 ,同时 数据 拟 合 线 也 更 贴近 于 1:1 线 。 


‘ 
有 
/| 


图 4 QYZ 四 组 模型 NEE 半 小 时 尺度 模拟 值 与 实测 值 的 RMSE 


3 结论 


本 文 以 千 烟 洲 亚热带 人 工 针 叶 林 和 长 白山 温带 针 阔 混交 林 生 态 系 统 为 例 , 应 用 模型 数据 融合 方法 为 不 同 
的 碳 循环 过 程 响应 模型 优化 参数 ,并 敌 选 最 适 模型 ,系统 分 析 了 水 分 对 不 同 生态 系统 碳 通 量 模拟 的 影响 。 主 
要 结论 如 下 : 

(1) 优 化 后 的 光合 .呼吸 参数 均 能 被 NEE 实测 数据 良好 约束 ,参数 后 验 为 正 态 分 布 标准 差 较 小 。 光 合 参 
数 的 生态 系统 差异 明显 :CBS 生长 季 的 LUE 和 4,, 均 高 于 QYZ, 但 参数 的 模型 间 差异 微小 ;呼吸 参数 的 生态 
系统 差异 明显 ;CBS 生长 季 的 Bn 和 010 均 高 于 QYZ, 同 时 参数 的 模型 间 差异 也 显著 :未 考虑 VPD 的 模型 高 估 
了 QYZ 的 Qi 值 (37%)、 低 估 了 QYZ 的 B 值 (45%)。 光 合 \ 呼 吸 参数 的 生态 系统 差异 准确 反映 了 温带 森林 与 
亚热带 森林 在 生长 季 的 碳 循环 特征 。 

(2) 同时 考虑 水 分 对 光合 .呼吸 作用 影响 的 模型 和 仅 考虑 水 分 对 光合 作用 影响 的 模型 分 别 是 QYZ 和 CBS 
的 最 优 模型 。 考 虑 水 分 因子 对 CBS 模型 的 模拟 性 能 并 无 明显 提升 ,对 QYZ 却 显著 提高 ,其 中 考虑 VPD 对 光 
合作 用 的 影响 是 QYZ 模拟 效果 提升 的 主要 原因 。 

(3) 不 同 的 水 分 响应 模型 在 CBS 碳 通 量 组 分 模拟 差异 较 小 ,而 在 QYZ 差异 明显 :未 考虑 水 分 的 模型 在 生 
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图 5 QYZ 干旱 年 生长 季 半 小 时 尺度 NEE 模拟 值 与 观测 值 散 点 图 
Fig.5 Half-hourly NEE from observations and the four models at QYZ 


长 季 高 估 了 GEP 总 量 2.0%(21.85 g C]m2? ) ,同时 更 大 幅度 地 高 估 了 RE 总 量 4.4%(38.02 g C]m2 ) ,从 而 导致 
整体 上 NEE 总 量 低估 于 实测 值 7.8% (18.55 g CAm ) 。 其 中 ,未 考虑 Sw 因子 主要 引起 RE 高 估 , 而 未 考虑 
VPD 不 仅 引起 GEP 高 估 , 也 导致 RE 高 估 。 同 时 ,在 QYZ 干旱 年 VPD 因子 对 碳 循环 的 影响 更 为 显著 ,未 考虑 
VPD 的 模型 不 能 模拟 出 干旱 期 NEE 日 变化 的 非 对 称 趋势 并 低估 了 NEE 的 峰值 。 

(4) 本 研究 中 提出 的 集成 MCMC 与 BIC 的 模型 数据 融合 方法 可 进一步 应 用 于 复杂 的 陆地 生态 系统 过 程 
模型 ,为 改进 生态 系统 模型 结构 .降低 模拟 的 不 确定 性 提供 了 有 效 途 径 。 
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